Antisite effect on ferromagnetism in (Ga,Mn)As 



O 
O 
(N 

0^ 



C/3 



R. C. Myers, 1 B. L. Sheu,^ A. W. Jackson, ^ A. C. Gossard,^ P. Schiffer,^ N. Samarth,^ and D. D. Awschalom^ 

^Center for Spintronics and Quantum Computation, 
University of California, Santa Barbara, CA 93106. 
^Department of Physics and Materials Research Institute, 
The Pennsylvania State University, University Park, PA 16802. 
(Dated: February 6, 2008) 

We study the Curie temperature and hole density of (Ga,Mn)As while systematically varying 
the As-antisite density. Hole compensation by As-antisites limits the Curie temperature and can 
completely quench long-range ferromagnetic order in the low doping regime of 1-2% Mn. Samples 
are grown by molecular beam epitaxy without substrate rotation in order to smoothly vary the As 
to Ga flux ratio across a single wafer. This technique allows for a systematic study of the effect of 
As stoichiometry on the structural, electronic, and magnetic properties of (Ga,Mn)As. For concen- 
trations less than 1.5% Mn, a strong deviation from Tc oc p"'^^ is observed. Our results emphasize 
that proper control of As-antisite compensation is critical for controlling the Curie temperatures in 
(Ga,Mn)As at the low doping limit. 

PACS numbers: 75.50.Pp, 75.70.-i, 71.55.Eq, 74.62.Dh 
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I. INTRODUCTION 

In (Ga,Mn)As, currently the most studied dilute mag- 
netic semiconductor (DMS), a long-range ferromagnetic 
interaction between dilute Mn^+ spins is mediated by va- 
lence band (or defect band) hole spinsiS*^ Most experi- 
mental work in this system has focused on achieving max- 
imal Curie temperatures (Tq's) in a technological drive 
to bring ferromagnetism in (Ga,Mn)As to more practical 
temperatures.^ In this vein, it was determined that a cru- 
cial limiting parameter of the Tq of as-grown (Ga,Mn)As 
is the incorporation of hole compensating crystalline de- 
fects, which are unintentionally formed due to the low 
growth temperature (~250 ° C) and low solubility of Mn 
in GaAs. One such defect, the Mn interstitial (Mui), is 
considered to be a double donor compensating two holes 
each^ These defects have shown a surprisingly high dif- 
fusivity allowing for their removal by using post-growth 
annealing at relatively low temperatures (~200 ° C) that 
drive them to the surface.^ Consequently, the three di- 
mensional hole density (p) and Tc after annealing can 
be drastically increasedfZ> Another well-studied defect in 
low temperature (LT) grown GaAs is excess arsenic. For 
substrate temperatures less than 300 ° C, excess arsenic 
is incorporated by occupying a Ga-site, called an antisite 
defect (AsGa) which - like the Muj defect is also a double 
donor The density of these defects can reach 10^" cm"'^, 
~1% of Ga-sites, which is the same order of magnitude as 
typical Mn doping levels j2. By post-growth annealing at 
temperatures above 500 ° C, Asca can be removed by the 
formation of metallic As precipitateSfiS. however at these 
temperatures Mn also precipitates out of the (Ga,Mn)As 
solid solution as MnAs nano-particles.^^ 

Most experimental work on (Ga,Mn)As has been 
aimed at maximizing the Tq by controlling the compen- 
sating defects in the high doping regime, i.e. Mn ^^5-9%. 
Such studies mainly focus on the reduction of Mni defects 
through post-growth annealing while loosely treating the 



problem of As non-stoichiometry. Asca defects are usu- 
ally ignored in the high doping regime since Muj defects 
are the main source of carrier compensation, and their 
removal has led to dramatic increases in p and Tc^ On 
the other hand, hole compensation due to Asca could 
be as large as 10^° cm~'^, which would shift Tq in this 
Mn-doping range by up to 10 K. 

In contrast, near the ferromagnetic threshold of ~ 1% 
Mn, compensation due to AsQa could dramatically al- 
ter the Tc's since here the defect and dopant densities 
are roughly equal. Additionally, the concentration of Mn 
necessary for the onset of ferromagnetism should be a 
strong function of carrier compensation since the ferro- 
magnetism depends strongly on p; in the low doping limit 
ferromagnetism is likely extrinsically limited by hole com- 
pensating defects. We also note that theoretical work pre- 
dicts a change in the mechanism of ferromagnetism in the 
low doping limit, near the insulator-metal transition^^iiS 
compared with the metallic limit at which Tc oc p'^-^^ as 
predicted by the Zener-model^ 

The effect of excess As on the lattice constant of 
(Ga,Mn)As was treated by Schott et al.,^^ who used a dif- 
ferent As-flux than another study^"* and noted a change 
in the extrapolated lattice constant of zinc blende MnAs 
due to AsQa incorporation. They foUowed-up their study 
by using 2 different As:Ga beam flux ratios, differing 
by 500% and noted large changes in the lattice con- 
stant of (Ga,Mn)As over a broad range of Mn-doping 
and substrate temperatures."^-^ These results were elabo- 
rated upon by Sadowski and Domagala, who measured 
lattice constant variation in Gao.96Mno.04As grown using 
3 different As:Ga flux ratios and noted an interplay be- 
tween Muj and AsQa incorporation for this Mn doping 
regime^ Campion et al. reported an improvement in 
structural and magnetic properties of (Ga,Mn)As films 
grown using As2 instead of As4.^^ Their paper, however, 
does not describe the fluxes used for each As-species nor 
how the effective III-V ratio varies between these cases. 
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The same group has reported impressively high Tc's over 
a broad range of Mn-doping levels, mentioning the im- 
portance of minimizing the As:Ga flux ratio although no 
data on the variation of this parameter were provided. 
Avrutin et al. carried out the first study that measured 
both the structural, electrical, and magnetic properties 
of (Ga,Mn)As using 2 different As:Ga flux ratios. They 
noted an optimization of the hole density, conductivity, 
and Tc for the low As:Ga flux ratio, where fewer Asoa 
should be incorporated. 

In this paper, we present the first systematic investi- 
gation of the effect of AsQa incorporation on the struc- 
tural, electronic, and magnetic properties of (Ga,Mn)As 
in the low Mn-doping regime. The As:Ga flux ratio 
is smoothly varied within individual samples by using 
molecular beam epitaxy (MBE) without substrate ro- 
tation. The geometry of the MBE chamber inherently 
provides a gradient in the As:Ga flux ratio across the 
substrate of ±50%, while maintaining a roughly uniform 
Mn flux. In order to avoid formation of Mni and to 
achieve a high sensitivity to carrier compensation, we 
chose Mn doping concentrations of 1 - 2%, in the re- 
gion of the onset of ferromagnetism. At Mn concentra- 
tions less than 1.5%, Muj has a larger formation energy 
than MuGa-"^ Thus no interstitial Mn will form and all 
hole compensation should originate from Asca- Film 
morphology is measured in- situ by reflection high en- 
ergy electron diffraction (RHEED) and ex-situ by atomic 
force microscopy ( AFM) . Out-of-plane epitaxial strain is 
measured by high-resolution x-ray diffraction (HRXRD). 
Electronic properties at room temperature are measured 
by Hall effect in the Van der Pauw and Hall-bar geom- 
etry, and magnetic properties are measured by SQUID 
magnetometry (Quantum Design MPMS). 

The epitaxial strain due to the incorporation of AsQa 
can be controlled by the As:Ga flux ratio, as previously 
reportediSii& More interestingly, the electronic and mag- 
netic properties reveal a strong sensitivity to Asca com- 
pensation with conductivities and hole densities (p) span- 
ning over two orders of magnitude leading to a transition 
from paramagnetism to ferromagnetism for constant Mn 
densities (Sec. IIIII and IVl)l . Hole compensation due to 
AsGa reduces p and therefore Tc- The As:Ga flux ra- 
tio is also observed to significantly alter film morphology 
(Sec. 13). In the As-rich condition, films attain max- 
imum smoothness, show a high degree of Asoa related 
carrier compensation, and exhibit reduced or suppressed 
Curie temperatures. For As:Ga flux ratios at which the 
p and therefore Tc are maximized (stoichiometric condi- 
tion, AsGa ~ 0), the films remain two-dimensional (2D) 
but are rougher than the As-rich material. In the Ga- 
rich condition (small As:Ga flux ratios), films become 
three-dimensional (3D) with the accumulation of excess 
Ga on the surface leading to an order of magnitude in- 
crease in film roughness. At the minimum As:Ga flux ra- 
tio investigated, the roughest surfaces are observed with 
hemispherical Ga-droplets. The group-Ill enriched sur- 
face suppresses Mn incorporation during growth, leading 



to lower p and Tq values than in the As- rich condition. In 
the low Mn-doped regime, where As-compensation limits 
hole density and ferromagnetic ordering, a strong devia- 
tion from the usual Tc oc p^-^^ behavior is observed with 
Tc oc p'^-^^ for X < 1.5%. Our results emphasize that 
proper control of hole compensation due to Asca is criti- 
cal to maximizing the value of Tc's in the low Mn-doping 
regime. The combinatorial MBE growth approach used 
here could be applied to more complex LT grown GaAs- 
based heterostructures in which the As:Ga flux ratio is a 
strong control parameter. 



II. MOLECULAR BEAM EPITAXY (MBE) 
GROWTH 

All samples are grown on 2-inch diameter GaAs (001) 
substrates in a Varian Gen-II MBE system manufac- 
tured by Veeco Instruments, Inc. with a source to sub- 
strate distance of 5.4-inches. The substrate tempera- 
ture is monitored by absorption band edge spectroscopy 
(ABES), which measures the substrate temperature in 
real-time through white-light transmission spectroscopy. 
The temperature stability during LT (Ga,Mn)As and 
GaAs growths is typically ±5 ° C . The growth rate of 
GaAs at the center of the wafer is 0.7 ML/sec as cal- 
ibrated by RHEED intensity oscillations of the specu- 
lar spot at 580 ° G. The Mn doping density is inferred 
from RHEED intensity oscillation measurementaSS of the 
MnAs growth rate performed at 240 ° C. Measurements 
of the Mn density by secondary ion mass spectroscopy 
(SIMS) indicate that the RHEED cahbration underesti- 
mates the density by up to 0.4% Mn, which provides a 
rough error estimate for the quoted Mn concentrations. 
We use As2 for all samples presented here, employing a 
valved cracker/sublimator. The beam equivalent pres- 
sure (BEP) of As2 and Ga are measured with a nude ion 
gauge at the center position of the substrate just prior 
to growth. The ratio of As2 and Ga BEPs, defined as 
As : Ga, is proportional to the ratio of the atomic fluxes 
of As and Ga delivered to the growing surface. Unfor- 
tunately, a reliable quantitative conversion from BEP to 
atomic flux is difficult due to the uncertainties in the ion- 
ization efficiencies of each beam speciesiSi A rough esti- 
mate gives that the atomic flux ratio is 2.5 times smaller 
than the BEP ratios, where we assume that only As2 are 
present with a temperature equal to that of the cracking 
zone. 

The following procedure is used for both non-rotated 
and rotated growths prior to deposition of the low tem- 
perature (LT) films. The wafer is first heated to 635 ° C 
with As : Ga ~ 35 for oxide desorption, monitored by 
RHEED, and subsequently cooled to 585 ° G. While ro- 
tating at 10 RPM, a 300-nm GaAs buffer layer is grown 
using 5-s growth interrupts for every 15 nm of deposi- 
tion resulting in a streaky 2x4 surface reconstruction 
RHEED pattern. While coohng the substrate to 250 ° G, 
the As- valve setting is reduced at 400 ° C for the desired 
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As : Ga value for LT growth. Prior to LT growth, a 
streaky 4x4 surface reconstruction is observed in the 
RHEED pattern. For rotated growths, once the sub- 
strate stabihzes at 250 ° C the LT growth begins. For 
non-rotated growths, the substrate rotation is stopped 
at 400 ° C and the wafer [110] crystal axis (perpendicu- 
lar to the wafer major flat) is aligned along the As : Ga 
gradient direction (y-axis). This can be reproducibly ac- 
complished noting the radial angle between the RHEED 
e-beam and substrate, fixed by the MBE chamber geom- 
etry, and finding a high symmetry crystal axis from the 
surface reconstruction pattern. All LT films are 100-nm 
thick at the center wafer position. After completion of 
the LT layer, they are cooled in a consistent manner to 
room temperature. 



III. ROTATED (GA,MN)AS: CARRIER 
COMPENSATION AND QUENCHED 
FERROMAGNETISM 

We first consider the case of rotated growth, examining 
the effect of Mn concentration and As : Ga on the elec- 
tronic and magnetic properties. Four samples are grown 
with the same Mn concentration (2.3%), but with various 
As : Ga. Hole density data are from samples measured 
in the Van der Pauw geometry in a magnetic field up 
to 0.2 T. Figure^ plots p and Tc of these samples as 
a function of As : Ga. Values of Tc are obtained from 
magnetization versus temperature scans as shown in Fig. 
^ selecting the point of maximum second derivative in 
M(T). The scans are performed during a field warm (50 
Oe) after a field cool (IT). The hole density decreases 
exponentially with the As : Ga, indicating that excess 
As in the form of Asca donors are compensating MuQa 
acceptors. Correspondingly, Tq decreases with p, and 
for As : Ga > 18 ferromagnetism is suppressed. This 
strong compensation effect is due to the fact that the 
AsGa donor and Mnca acceptor densities are similar in 
magnitude in this low Mn-doping regime, Mn < 2 %. 

Using As : Ga = 10.3, for which the largest p and 
Tc are obtained (Fig. ^), a series of four additional 
samples are grown with varying Mn concentration. For 
this series, variation in p and Tc occurs due to the varying 
concentration of Mnca, Fig. ^p. Clearly the hole density 
does not increase linearly with Mn concentration, which 
is indicative of a carrier compensation threshold. Based 
on this observation, a more careful optimization of As : 
Ga is necessary, as discussed in the next sections. 

Figure ^ plots Tc as a function of p, summarizing 
data from Fig. ^ a and b. We observe that varying the 
concentration of compensating Asca donors (by varying 
As : Ga) has a qualitatively identical effect as varying the 
concentration of magnetic MuQa acceptors with respect 
to the carrier density and magnetic transition tempera- 
ture. A fit of Tc oc p°-33 to the ferromagnetic samples in 
Fig. nj;, Mn doping > 1.5%, reveals the usual functional 
dependence/^ 
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FIG. 1: Quenched ferromagnetism by As-compensation. 
Data from (Ga,Mn)As films grown at 250 ° C using substrate 
rotation, (a) Samples were grown using a constant growth 
rate and Mn concentration but varying As : Ga. Room tem- 
perature hole density (p) is plotted on the left vertical axis as 
a function of As : Ga. The corresponding Curie temperature 
(Tc) is plotted on the right vertical axis, (b) Similar data 
plotted for samples grown with a constant As : Ga, but vary- 
ing Mn concentration, (c) Data in (a) and (b) are replotted 
with Tc as a function of p. A fit to Tc shown for 

the ferromagnetic samples, (d) Magnetization versus temper- 
ature scans (field-warmed with B = 50 Oe) used to extract 
Tc shown in (b). 
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IV. NON-ROTATED (GA,MN)AS: 
COMBINATORIAL VARIATION OF As 



Ga 



In order to systematically study the effect of the ex- 
cess As on the electronic and magnetic properties of 
(Ga,Mn)As, we use a method of non-rotated growth in 
which the geometry of the MBE system provides a con- 
tinuous variation in As : Ga across the wafer^ Thus 
the advantage of a combinatorial approach is employed 
in which a single non-rotated growth provides the equiv- 
alent of twenty rotated growths of various As : Ga. The 
phase space of As : Ga and Mn concentration can be 
mapped-out in high-resolution while performing a much 
smaller number of growths than otherwise required. 

Figure |2t schematically shows the position of the Ga 
Knudsen cell and the As valved-cracker/sublimator rela- 
tive to the substrate position. The angle between the Ga 
and As2 molecular beams provides a continuous variation 
in As : Ga across the y-axis of the substrate, as defined 
in the schematic. The variation of the As flux is mea- 
sured by depositing a thick As film on a quartz wafer, 
covered by a shadow mask along the y-axis, with the 
substrate stabilized below room temperature. We have 
assumed that the deposition rate of As under these con- 
ditions is proportional to the As flux. A profilometer is 
used to measure the thickness of the As film along the y- 
axis. The normalized flux variation of As, plotted in Fig. 
Eb, is calculated by normalizing the thickness along y by 
the thickness at y = 0. The variation of the Ga-flux is 
measured using an optical measurement of a distributed 
Bragg reflector (DBR) and cavity structure grown in 2 
separate steps. First, a 10 period GaAs/AlGaAs DBR is 
grown on a rotated 2-inch substrate in a larger, more uni- 
form system (Veeco Gen III) with a source to substrate 
distance of 11-inches. This results in an extremely uni- 
form DBR which is As-capped and transferred into the 
Gen II system. The As-cap is desorbed and a GaAs opti- 
cal cavity layer is grown on the sample without substrate 
rotation. The reflectance spectrum of the sample is mea- 
sured ex-situ at different points across the wafer and is 
fit to an optical model of the reflectance. Because the 10 
period DBR was almost completely uniform, the thick- 
ness of the non-uniform GaAs layer grown in the Gen 
II could be accurately and unambiguously determined 
by fitting the measured optical spectrum to a computer 
model of the structure. Since the deposition is performed 
at 585 ° C, the growth rate is independent of the As over- 
pressure, being limited by the group-Ill flux only. The 
normalized Ga rate, proportional to the Ga-flux, along 
the y-axis is plotted in Fig. I^J- Both the As and Ga flux 
variation are well flt by a second order polynomial (lines), 
which are used to convert y-axis position for non-rotated 
growths to As : Ga. 

Next we consider the Mn flux variation across the sub- 
strate. Since the Mn Knudsen cell is roughly perpendicu- 
lar to the y-axis, variation along this direction is expected 
to be small. In order to measure the Mn variation across 
the substrate, a l-/zm thick GaAs:Mn layer with Mn ~ 
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FIG. 2: Geometric flux gradients in the MBE system. 
As : Ga is continuously varied and optimized in a single non- 
rotated growth, (a) Deposition geometry (not to scale) in 
the MBE system. The relative positions of the As, Ga, and 
Mn Knudsen cells with respect to the substrate are shown 
from side and front views, (b) The flux gradient of Ga and 
As molecular beams along the y-axis of the wafer normalized 
to the center flux value. Measured values (data points) and 
flts to second order polynomials (lines) are shown, (c)-(d) The 
Mn:Ga flux variation along the x-axis (c) and y-axis (d) of the 
wafer measured from hole density (p) variation in non-rotated 
samples grown at 400 C, where As-flux does not affect carrier 
density. The line in (d) is the inverse of the Ga-flux variation, 
shown as a line in (b), which accounts for the variation in the 
Mn:Ga ratio along the y-axis. 



4x10^'' cm~^ is deposited at 400 ° C with an As2 over- 
pressure. At these conditions Mn incorporation is fully 
substitutional and excess As incorporation is negligible 
(< 1 X 10^^ cm~^), thus the hole density is proportional 
to the Mn:Ga flux ratioi^^ The hole density variation 
across the x and y axes due to Mn:Ga variation is plot- 
ted in Fig. [21c and d, respectively. Along either axis, the 
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Mn:Ga ratio varies by 5-10%, which would correspond to 
variation in the Mn concentrations of 0.05 to 0.1%. As 
will be clear (Sec. IVIL this variation in Mn concentration 
is negligible with respect to the As compensation effect 
under investigation. 



V. FILM MORPHOLOGY AND STRAIN 
VERSUS As : Ga 

To study the effect of As : Ga on the film morphology, 
the RMS roughness of non-rotated samples is measured 
across the y-axis by AFM. We describe results from two 
non-rotated films, without Mn doping and with 1.5% Mn. 
Both 5 X 5-/im^ and 1 x l-/im^ area scans are measured ev- 
ery 3 to 4 mm along the y-axis of both wafers. The RMS 
roughness from these scans is plotted as a function of 
As : Ga (converted from y-position from Fig. calibra- 
tion) and representative AFM images of each region are 
shown in Fig. U] 

In the As-rich region from 11 < As : Ga < 14, AFM 
images reveal fully coalesced films with RMS roughness of 
0.2 nm or less indicative of a 2D growth mode. RHEED 
observations of such films, from rotated growths at the 
same As : Ga condition, show sharp and streaky 2D re- 
construction. In the As-rich to stoichiometric region from 
9 < As : Ga < 11, films show a nanoscale- granular struc- 
ture. For decreasing As : Ga, this granular structure 
becomes more distinct and the film roughness increases, 
however it remains below 1 nm RMS. RHEED observa- 
tions for films grown in this As : Ga range show a spotty 
reconstruction pattern indicative of a rough, but 2D sur- 
face, which we have found to be characteristic of films 
near the stoichiometric condition. In this region, the hole 
density of Mn-doped films reaches a maximum indicating 
the minimum of compensating Asca defects (Sec. IVIjl . 

Just below As : Ga ^ 9, the RMS film roughness in- 
creases by a factor of ten as As : Ga transitions from 
the stoichiometric condition to Ga-rich. AFM images 
of this region show the formation of 10-nm tall clover- 
shaped features 60-nm in diameter with a crater in the 
center. The RHEED pattern of these films shows addi- 
tional spots indicative of 3D features on the surface. For 
As : Ga < 9, the surface roughness continues to increase 
with the appearance of larger droplet-like features. At 
the lowest As : Ga ratio, these hemispherical droplets are 
50-nm tall and 240-nm in diameter. The large change in 
roughness across the non-rotated wafers is observable by 
eye under intense illumination, where the surface transi- 
tions from mirror finish at the As-rich region to a hazy 
surface in the Ga-rich region. 

As observed in the roughness plot of Fig. El it appears 
that Mn-doping suppresses the onset of roughening to- 
ward smaller As : Ga. In contrast, since Mn-doping in- 
creases the effective group-Ill flux, one would expect a 
shift in the onset of roughness to higher As : Ga. This 
smoothing may be due to a surfactant effect of Mn on the 
surface. As discussed in Sec. I VII the highest hole densi- 
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FIG. 3: Morphology of GaAs and (Ga,Mn)As films as a func- 
tion of As : Ga (non-rotated growths). AFM micrographs 
(1x1 fira^ scans) are plotted for a GaAs growth, but are qual- 
itatively identical to (Ga,Mn)As films in the same wafer po- 
sition. The RMS roughness measured from such AFM scans 
is plotted as a function of As : Ga (along y-axis of wafer) 
for both a GaAs and (Ga,Mn)As sample. The stoichiometric 
region is shaded grey, where excess As is minimized, two- 
dimensional growth (2D) is maintained, however, some gran- 
ular structure is observed. 



ties are observed in (Ga,Mn)As grown with As : Ga just 
before the transition from 2D to 3D surface, indicating 
that the growth kinetics at this stoichiometric condition 
do not significantly increase the density of hole compen- 
sating Mui donors. 

HRXRD scans are obtained using a triple-axis (Phillips 
X'Pert MRD Pro) thin film diffractometer. We mea- 
sure uj-29 scans near the (004) Bragg peak of the (001) 
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FIG. 4: Film strain due to As non-stoichiometry and Mn in- 
corporation in GaAs and (Ga,Mn)As non-rotated films mea- 
sured by u-26 HRXRD scans near the (004) peak, (a)-(b) 
Data from a GaAs film at the As-rich wafer position (a) and 
stoichiometric position (b). Data (points) and dynamical sim- 
ulations (lines) are plotted for various Asoa densities, (c)-(d) 
Data from a (Ga,Mn) As film at the As-rich wafer position (c) 
and stoichiometric position (d). Data (points) and dynamical 
simulations (lines) are plotted using different Mn densities. 



GaAs substrate. The diffraction patterns obtained for 
the non-rotated LT GaAs film at the position of highest 
As-compensation {As : Ga = 13 : 1) and at the stoi- 
chiometric condition (As : Ga = 9:1) are plotted in 
Fig. ^a and b, respectively. In the As-rich position, the 
epilayer peak appears as a shoulder at lower angle, re- 
flecting the increase in the out-of- plane lattice constant, 
due to AsQa incorporation. We also note the thickness 
fringes indicative of the high quality of the film surface 
and interface. A dynamical HRXRD program is used to 
simulate the data™ assuming that the lattice constant of 
(Gai_z,Asz)As, a = 5.654 -I- 1.356z.^ Best fit simulations 
in Fig. 2t give a film thickness 90 nm and Asca = 4x10^^ 
cm^'^. The additional fit with Asca = 5x10^^ cm~'^ 



shows the error in the Asca estimate is about 1x10"'^^ 
cm~'^. The thickness of the epilayer is 10-15% thinner at 
the top edge than at the center of the wafer due to the 
smaller Ga-fiux at this position, as shown in Fig. ^p. At 
the stoichiometric position the epilayer shoulder becomes 
small, only visible as a slight asymmetry on the left side 
of the substrate diffraction peak in Fig. ^p. Simulations 
using a film thickness of 100 nm yield Asca = 4x10^^ 
cm~'^, which is smaller than the uncertainty of the fits. 
A more precise estimate of the Asca density is found 
from the hole density variation, discussed in Sec. IVII 
Thickness fringes disappear since the lattice constant of 
the epilayer is almost indistinguishable from that of the 
substrate. 

Figure 01 c and d plot the diffraction patterns from a 
non-rotated (Ga,Mn)As film with 1.5% Mn from the As- 
rich and stoichiometric positions, respectively. For these 
layers, strain due to Mn incorporation dominates the epi- 
layer lattice constant obscuring the Asca effect. Simula- 
tions are performed assuming that the lattice constant 
of (Gai_^,Mn^)As, a = 5.654 + 0.476x.i3 The fits in- 
dicate a film thickness of 85 nm and 95 nm, and Mn 
concentration of 1.5% and 1.3% at the As-rich and sto- 
ichiometric positions, respectively. The variation in Mn 
concentration is not explained by the Mn:Ga ratio varia- 
tion along y. Fig. which would predict a decrease in 
Mn concentration of up to 0.1%. The additional strain 
reduction may be due to the minimization of Asca de- 
fects at the stoichiometric position, but the uncertainty 
of the HRXRD simulations prevents a more quantitative 
analysis. 



VI. CARRIER DENSITY AND CURIE 
TEMPERATURE VERSUS As : Ga 

A. Hole density gradient 

Figure \^ a and b show the variation in hole density 
for a series of non-rotated (Ga,Mn)As films of various 
Mn concentrations (a) along the y-axis, and (b) with y 
converted to As : Ga, using the Fig. |2b calibration. 
All the hole density data in the figures are from room 
temperature measurements of samples prepared in the 
Van der Pauw geometry (3x3 mm^) in a magnetic field 
up to 0.2 T. Selective measurements of the hole density 
of patterned Hall-bars in magnetic fields up to 2 T agree 
with the Van der Pauw measurements to within 50%. 
Additional high field measurements at 7-9 T deviate from 
the low field data by up to 50 %. 

The hole density varies by over two orders of magni- 
tude across each wafer showing a maximum value near 
y ^ 6 mm corresponding to As : Ga ~ 9. The results 
are explained by a decreasing number of AsQa donors 
that compensate MuQa acceptors as As : Ga is reduced 
from the As-rich condition. The maximum of p along 
y corresponds to the position at which Asca are mini- 
mized, the stoichiometric condition. A further reduction 
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FIG. 5: The As-compensation gradient in non-rotated 
(Ga,Mn)As films. The room temperature hole density (p) 
variation along films grown with different Mn concentrations, 
(a) as a function of y-position, and (b) as a function of 
As : Ga. (c)-(d) Similar data for samples with the same 
Mn concentrations but different Aa : Ga at y = 0. Lines in 
each plot connect data points from the same wafer. 



in As : Ga leads to a Ga-rich surface that limits Mn in- 
corporation, and leads to a ten-fold increase in surface 
roughness marking a transition from 2D to 3D film mor- 
phology (Sec. 0). The y-position of the maximum in 
p matches in four separate non-rotated growths demon- 
strating the reproducibility of this technique. When the 
Mn density is reduced, the hole-density peak narrows on 
the As-rich side, indicating that samples with lower Mn- 
doping density are more sensitive to the AsQa compensa- 
tion effect. This follows if As : Ga controls the density of 
hole compensating defects independently from the dop- 
ing density of Mn, as expected for this doping regime. 

We note that for the 0.75% Mn doped wafer, when 
As : Ga ^10 the hole compensation is nearly complete, 
p ~ 10^^ cm~^. For As : Ga ^ 9, the compensation is 
minimized, p ~ 10^^ cm~'^. In this case, a 10% reduc- 
tion in As : Ga reduces the Asoa compensation by three 
orders of magnitude. This demonstrates the strong sen- 
sitivity of the AsGa incorporation to As : Ga, and the 
high degree of As-flux control, within 10%, necessary to 
properly minimize this defect. Previous studies of the 



effect of As : Ga ratio reported only two or three values 
of this parameter in increments of 200% or morei^'i^^ 

In the Ga-rich position {As : Ga < 9), the hole density 
decreases along with As : Ga, but remains constant as 
the Mn doping density is varied. Fig. [Sja. This behavior is 
explained by a solubility limit for Mn, which is below any 
of the Mn doping levels used here. Decreasing As : Ga 
below the stoichiometric value leads to a Ga-rich surface 
that serves as a kinetic barrier for the incorporation of 
group-Ill substitutional dopants, such as Mn. One might 
consider the alternative possibility that interstitial incor- 
poration of Mn might occur in these circumstances and 
could explain the reduction in hole density in the Ga-rich 
region. However, if Muj donors were incorporating and 
compensating the Mnca acceptors, then the carrier den- 
sity should change as the Mn-doping density is altered, 
which is not the case here. Preliminary SIMS data show 
a decreasing concentration of Mn in the Ga-rich region 
as the thickness increases, whereas the As-rich region re- 
veals a relatively constant Mn depth profiles throughout 
the Mn-doped layer. In the Ga-rich region the Mn den- 
sity can be ~ 10^" cm~^ lower than in the As-rich region. 

The position of the hole density peak can be controlled 
by adjusting As : Ga at y = 0, which shifts the position 
of the stoichiometric condition along y. To demonstrate 
this, two wafers with identical Mn concentration of 1.5 
% are grown but with different As : Ga at y = 0. The 
variation of p along y is plotted in Fig. By increasing 
As : Ga at y = from 8:1 to 9:1, the peak position (sto- 
ichiometric condition) is shifted by ~ 1 cm toward the 
center of the wafer in a direction away from the As source. 
This shift is observable by eye as a change in the posi- 
tion of the transition region, where the surface changes 
from mirror finish to hazy. Converting the y-position to 
As : Ga should fully overlap the hole density peaks. Fig. 
131. The peaks do not completely overlap indicating that 
additional factors must be considered, such as the varia- 
tion in the growth rate across non-rotated wafers, which 
may play a role in the density of compensating Asca de- 
fects. In the present analysis we have considered only the 
effect of As:Ga. 

Finally, we compare samples of equal Mn doping den- 
sity grown with and without rotation. The hole densities 
of the rotated samples. Fig. match the hole densities 
for the non-rotated samples with similar Mn doping in 
the As-rich condition of As : Ga ~ 10.5, Fig. |5|d. This 
indicates that the rotated sample set, presented in Fig. 
^P, is grown under a slightly As-rich condition, on the 
right side of the hole density peak of Fig. ^d. 



B. Curie temperature gradient 

The Curie temperature variation with As : Ga for two 
non-rotated (Ga,Mn)As films is plotted in Fig. (Sja and 
b, together with the hole density variation. Tq is deter- 
mined from scans similar to those plotted in Fig. ^i. Tq 
tracks the p variation with As : Ga, displaying a maxi- 
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FIG. 6: The Curie temperature gradient in non-rotated 
(Ga,Mn)As films due to As-compensation. The hole density 
(p) and Tc are plotted as a function of As : Ga for two wafers 
with different Mn concentrations, (a) 1.25%, and (b) 1.50%. 
Lines guide the eye. The stoichiometric region is shaded grey, 
where excess As is minimized. 



mum in Curie temperature in the same region where the 
hole density peaks. In the As-rich portion of the wafer, 
AsQa donors Umit the hole density and therefore Tq. This 
effect is very clear, for example, with As : Ga > 12, 
where ferromagnetism is completely suppressed in both 
1.5% and 1.25% samples due to hole compensation from 
As non-stoichiometry in the form of Asca donors. 

The width and shape of the Tc and hole density peaks 
are quite different. In general, the Tq peaks are wider 
than the hole density peaks on the As-rich side, but are 
sharper than the hole density peak on the Ga-rich side of 
the wafer. To explain this behavior, we replot data from 
several of the non-rotated (Ga,Mn)As wafers to examine 
the dependence of Tc on the hole density. This depen- 
dence is plotted for the As-rich region {As : Ga > 9) in 
Fig. [7^, and for the Ga-rich region {As : Ga < 9) in Fig. 

In the As-rich region, ferromagnetism is observed over 
a broad range in hole density from 10^^ to 10^^ cm^'^. 
For Ga-rich material, however, ferromagnetism is not 
observed over the same large range in hole densities 
where ferromagnetism is observed in the As-rich mate- 
rial, shaded grey. As described in Sec. Ivl and IVI Al hole 
density and film morphology measurements indicate that 
in the Ga-rich region Mn solubility is limited to well- 
below the level required for ferromagnetism, < 0.5 %. 
This is probably due to the kinetic disadvantage of a 
group-Ill dopant on a heavily Ga-rich surface because of 



FIG. 7: Curie temperature (Tc) versus hole density (p) in 
several non-rotated (Ga,Mn)As films showing deviation from 
p'''^^ . (a) Data are plotted from the As-rich region of the wafer 
{As : Ga > 9), and (b) from the Ga-rich region {As : Ga < 9). 
The shaded region marks p- values at which ferromagnetism is 
observed in As-rich material (a) and not in Ga-rich material 
(b). The line fit of data in (a) to Tc oc p"'"^ is replotted in 
(b). 



out-competition. The discrepancy in Fig. [7| between the 
Ga-rich region and As-rich region offers further support 
for this hypothesis. In this interpretation, Ga-rich condi- 
tions limit MuQa to low doping densities of less than 10^^ 
cm^'^ ( < 0.05% Mn). Whereas for similar hole densities 
in the As-rich material Mnca ~ 10^° cm"^ (--1% Mn). 
Thus ferromagnetism is possible in the As-rich material 
since there are sufficient Mnca acceptors, but the Tc is 
limited by hole compensation due to Asca donors. In 
the Ga-rich material, Asoa compensating donors are not 
present, but the Mnca density is well below the ferro- 
magnetic limit. 

The data in Fig. [TJi for which Tc ^ are fit to 
Tc oc (line). This fit describes the data from two 

non-rotated wafers with Mn = 1.5% (Fig. |S|;) and one 
wafer with 1.25% Mn over a large range in p and Tc, over 
13 data points. In the rotated growths presented in Fig. 
n pO.33 behavior is observed for samples with > 1.5% Mn. 
There is no discrepancy between the non-rotated and ro- 
tated (Ga,Mn)As films if we conjecture that a change in 
the exponent from 0.33 to 0.09 occurs as the doping is re- 
duced below 1.5% Mn. Since these doping levels are close 
to the ferromagnetic limit, a change in the dependence 
of ferromagnetism is not unexpected. The impurity band 
is fully merged with the valence band in the metallic re- 
gion of high Mn-doping, where the pO -^^ jg observed. But 
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at low Mn-doping levels close to the insulating regime, 
the impurity band is separated from the valence band 
(Mott gap) where a disorder-based model of ferromag- 
netism is probably more accurate. Erwin and Petukhov 
predicted that a strong deviation fr'om the metallic be- 
havior of Tc cx should occur in the low Mn-doping 
regime, ^^1.5% Mn, with Tc reaching a maximum in a 
partially compensated condition (lower p) and ferromag- 
netism disappearing for zero compensation (highest p); 
their model assumed that compensation was due to Muj 
alone^^ This behavior is not clearly seen in our data, how- 
ever the observed AsQa compensation was not treated in 
the model of Erwin and Petukhov. 



VII. CONCLUSIONS 

The electronic and magnetic properties of (Ga,Mn)As 
reveal a strong sensitivity to Asca donors that compen- 
sate Mnca acceptors. By varying As : Ga, the hole den- 
sities of (Ga,Mn)As films are varied by more than two 
orders of magnitude, and, as a result, fcrromagnetism 
can be enhanced or completely suppressed. For As-rich 
growth, films show the maximum smoothness, but low- 
est p and Tc- If As : Ga is reduced, the film rough- 
ness increases up to the stoichiometric condition, where 
p and Tc are maximized. For Ga-rich growth, the film 
roughness increases by an order of magnitude as excess 
Ga accumulates on the surface. The group-Ill rich sur- 
face suppresses substitutional Mn incorporation leading 



to lower p and Tc values than in the As- rich condition. In 
the low Mn-doped regime, where As-compensation limits 
hole density and ferromagnetic ordering, a strong devi- 
ation from the usual Tc oc p^-^^ behavior is observed 
with Tc oc p'^ o^ for x < 1.5%. This dependence applies 
equally well to stoichiometric (Ga,Mn)As as it does to 
heavily Asca compensated (Ga,Mn)As. This change in 
the hole density dependence of Tc occurs in the doping 
region where the effects of a Mott gap and disorder are 
expected to alter the mechanism of ferromagnetismi^ii^ 
Our results emphasize that precise control of hole com- 
pensation due to AsQa is critical to maximizing the value 
of Tc's in the low Mn-doping regime. The combinato- 
rial MBE growth approach used here could be applied 
to more complex LT grown GaAs-based heterostructures 
in which As : Ga is a strong control parameter. The 
rapid and high-resolution mapping of growth parameter 
phase-space by this technique allows for systematic inves- 
tigations not practicable using more traditional methods. 
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